Finite-size effects on the dynamic susceptibility of CoPhOMe 
single-chain molecular magnets in presence of a static magnetic field 
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The static and dynamic properties of the single-chain molecular magnet [Co(hfac)2NITPhOMe] 
are investigated in the framework of the Ising model with Glauber dynamics, in order to take into 
account both the effect of an applied magnetic field and a finite size of the chains. For static fields of 
moderate intensity and short chain lengths, the approximation of a mono-exponential decay of the 
magnetization fluctuations is found to be valid at low temperatures; for strong fields and long chains, 
a multi-exponential decay should rather be assumed. The effect of an oscillating magnetic field, with 
intensity much smaller than that of the static one, is included in the theory in order to obtain the 
dynamic susceptibility xi^)- We find that, for an open chain with JV spins, x( w ) can be written as 
a weighted sum of N frequency contributions, with a sum rule relating the frequency weights to the 
static susceptibility of the chain. Very good agreement is found between the theoretical dynamic 
susceptibility and the ac susceptibility measured in moderate static fields (-ffdc < 2 kOe), where 
the approximation of a single dominating frequency turns out to be valid. For static fields in this 
range, new data for the relaxation time, r versus Hdc, of the magnetization of CoPhOMe at low 
temperature are also well reproduced by theory, provided that finite-size effects are included. 



PACS numbers: 75.50.Xx,75.40.Gb,75.10.Hk,76.20.+q 

I. INTRODUCTION 

R is commonly admitted that the plain Ising Hamil- 
tonian does not contain any dynamics^ In fact, when 
considering a system of Ising spins, <7i, localized at the 
sites, i, of a lattice 

n = -j I Y u °i°h tTi = ±i (i) 

the physically interesting quantities, <7j's, commute with 
rl. However, for a system in contact with a heat bath, 
a stochastic dynamics can be introduced by means of 
a master equation which assumes Markovian processes 
inducing random flips between different states. One of 
the few cases in which the problem can be solved an- 
alytically is for a one-dimensional lattice, zero external 
magnetic field, and an opportune choice of the transition 
probability, as devised by Glauber— some decades ago. 
He calculated the dynamic susceptibility within a linear 
response framework, and found that the uniform mag- 
netization decays exponentially, with a relaxation time 
given by an Arrhenius law 

T(T) = r e^, (2) 

where 1/tq is the relaxation rate for an isolated spin. 
Considering that, at low temperatures, the correlation 
length for model ([T]) in one dimension is given by £ oc 

2 J I 

e k B T * the dynamic critical exponent results in z = 2. 
Glauber's dynamics has been ever since applied in very 



different areas, comprising structural phase transitions-, 
neural networks^, chemical reactivity^ and even bio- 
socio-econo-phy sics 

The recent discovery of single-chain magnets (SCMs) 10 
spurred renewed interest for Glauber's dynamics in mag- 
netic nanomaterials. Such systems show magnetic hys- 
teresis without the onset of three-dimensional magnetic 
ordering. At very low temperatures, the relaxation of 
the magnetization is so slow that also other very inter- 
esting dynamic phenomena have been observed with un- 
precedented clarity, including collective reversal^ and 
quantum tunneling 1 ^. As a result, after the discovery 
of the archetypal SCM [Co(hfac)2NITPhOMe] (hereafter 
denoted CoPhOMe)^ the number of SCM compounds 
has been rapidly increasing* 1 ^ - — 

The strong exchange interaction and one-dimensional 
character of CoPhOMe 10 make it the ideal system where 
to observe the long-predicted 2 - slow relaxation of the mag- 
netization. In CoPhOMe, owing to the very high value of 
the exchange constant (J///cs ~ 80 K), the correlation 
length £ is huge, at low temperatures, in a zero magnetic 
field. Consequently, the unavoidable presence of even a 
small density of defects causes the chain to break into 
finite segments whose average length, L, can be much 
smaller than the correlation length, £. In such a finite- 
size regime, L <C £, the dynamics of the system in zero 
field is modifie d 21 1 22 with respect to Glauber's analysis of 
the infinite chain.— In particular the relaxation time, mea- 
sured in CoPhOMe using ac susceptibility and SQUID 
magnetization decay techniques, was found to follow an 
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Arrhenius law with a halved energy barrier— — 
t(T)=t (L) e^r, 



(3) 



in agreement with theoretical predictions of finite-size 
effects,^ 

These effects were systematically investigated in 
CoPhOMe by introducing non-magnetic impurities. 26 
For nominally pure and impure samples, the complex sus- 
ceptibility x( w ) — x'fa) + *x"( w ) was measured in pres- 
ence of a moderate static magnetic field (iJdc = 2 kOe) 
and of a much smaller ac field oscillating at frequency uj. 
A two-peaked structure was found in x'(w) as a function 
of temperature: the low temperature peak is frequency- 
dependent, while the high temperature one is not. On 
the basis of transfer matrix calculations for the static sus- 
ceptibility of the doped chain, the low temperature peak 
was attributed to finite-size effects^ Anyway, the fre- 
quency dependence of the peak for x'(cj) and x"(w) ^ 
remains unexplained. Static and dynamic susceptibilities 
with a similar behavior have also been observed in differ- 
ent cobalt-organic single-chain magnets,— in presence of 
a moderate static magnetic field {H& c — 0.5 kOe) and of 
lattice imperfections. The explanation of such features is 
thus becoming more pressing, as it can constitute an im- 
portant tool for the analysis of the properties of a whole 
class of magnetic systems. 

In systems with a very large correlation length, like 
the one-dimensional Ising model at low temperatures, 
the introduction of an external magnetic field can have 
dramatic consequences. A static field H& c strongly de- 
presses the correlation length £j£ and this fact, in turn, 
should also strongly affect the dynamic susceptibility. 
The Glauber dynamics of the infinite Ising chain model 
in an external magnetic field was studied some years ago 
in order to describe the kinetics of the helix-coil transi- 
tion in biopolymerS: 2 - For single-chain magnets, a the- 
oretical and experimental study was recently performed 
by Coulon et a/.— focusing on the relaxation time of the 
magnetization fluctuations. In addition to a static mag- 
netic field, their theoretical analysis 2 ^ considered finite- 
size effects, relevant in SCMs. A local-equilibrium ap- 
proximation was adopted 2 ^ in order to truncate the in- 
finite hierarchy of kinetic equations for finite open Ising 
chains in H^c ^ 0. The main advantages of this approxi- 
mation, first proposed by Huang 30 for infinite chains, are: 
(i) it provides the exact steady-state solution^ in con- 
trast with the mean field approximation; (ii) it is valid 
for any value of the applied field, in contrast with per- 
turbation methods^ 

In this work, we develop the theoretical framework nec- 
essary to analyze the ac susceptibility measurements of 
single-chain magnets in presence of a static magnetic field 
Hdc- We include in the theory 2 ^ the effect of an oscil- 
lating magnetic field with intensity much smaller than 
that of the static one, using a linear response framework. 
With these theoretical tools, we can account for the dy- 
namic behavior of CoPhOMe and other single-chain mag- 
nets. We directly compare the calculated behavior to 



previous 2 ^ and new experimental measurements of both 
the relaxation time and the ac susceptibility in presence 
of a static magnetic field i?dc- In this way, we can repro- 
duce the temperature and frequency dependence of the 
ac susceptibility of SCMs, measured in moderate static 
fields (-f/dc < 2 kOe) and in presence of crystal defects 
and/or non-magnetic impurities. For static fields in this 
range, new data for the relaxation time of the magneti- 
zation of CoPhOMe, r versus H^c, are analyzed at low 
temperature. They are well reproduced by theory, pro- 
vided that finite-size effects are included. 

The paper is organized as follows. In Section |n] we 
present both the real system and the simplified model 
that we adopt to catch the essentials of its stochastic dy- 
namics. In Section IlIII we calculate, for pure and doped 
chain systems, the temperature dependence of the mag- 
netization and static susceptibility in presence of a static 
magnetic field. In Section ITVl the theoretical framework 
for the calculation of the dynamic susceptibility in pres- 
ence of a static magnetic field is first developed for the in- 
finite chain, and then generalized to an open, finite chain 
with N spins; next we present and discuss the approxima- 
tion of a single dominating contribution, with character- 
istic frequency f2 c , to the dynamic susceptibility. In Sec- 
tion[V] the results of our explicit calculations, performed 
using parameters suitable for describing CoPhOMe, are 
shown and discussed. Finally, the conclusions are drawn 
in Section IVI1 and some technical details are reported in 
the Appendices, for the reader's convenience. 



II. THE MODEL 

The magnetic properties of CoPhOMe are determined 
by Co(II) ions, with an Ising character and effective 
S = 1/2, and by NITPhOMe organic radicals, magnet- 
ically isotropic and with s — l/2 i 10 ' 23 ' 31 The primitive 
magnetic cell is made up of three cobalt ions and three 
radicals. The spins are arranged in a helical structure, 
with the helix axis coincident with the crystallographic 
c axis of the chain. Although the effective spins of the 
two types of magnetic centers have the same value (1/2), 
the gyromagnetic factors are different. The gyromag- 
netic factor gn of the organic radical is isotropic, while 
cobalt is strongly anisotropic and gco = 3c ^ <?co' with 
<7Co 3> 9r- The Hamiltonian currently used^ 3 - to describe 
a CoPhOMe chain of N spins reads 



N/6 3 (- 

^ = - ^2 ^2 \ J ai - 2m 

1=1 m=l I 



C7,2m-1 + a l,2m+l 



fi B H 



(Z2m • e H ) >, (4) 



where a is the spin variable, I is the magnetic cell index 
and m the site label. The nearest neighbor exchange is 
antifcrromagnetic and rather strong, \ J\/ks ~ 80 K; en- 
denotes the direction of the applied magnetic field H; 
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finally, all local axes along which the spins are aligned, 
Zfe (k = 1, • • • , 6), form the same angle 6 « 55° with the 
c axis. It follows that the sublattice magnetizations are 
not compensated along c, whereas they are compensated 
within the plane perpendicular to the chain axis. Thanks 
to its ferrimagnetic and quasi-one-dimensional character 
(the ratio between interchain and intrachain exchange 
constants is less than 10 -6 )^ CoPhOMe was the first 
magnetic molecular compound to display slow relaxation 
of the magnetization at low temperatures for H = 0^2- a 
feature which was predicted long time ago by Glauber 2 
in a one-dimensional model of Ising spins, coupled by a 
nearest neighbor ferromagnetic exchange and interacting 
with a heat reservoir, causing them to change their states 
randomly with time. 

In this paper we are primarily concerned with analyz- 
ing finite-size effects on the spin dynamics of CoPhOMe 
in presence of a non-negligible static magnetic held. Since 
model Q is too involved, in the following we adopt a sim- 
plified model, yet able to catch the essentials of the dy- 
namic behavior. Namely, we make the approximation of 
an open Ising chain with N equal spins, all with a = ±1 
and the same (isotropic) gyromagnetic factor g, coupled 
by an effective nearest neighbor ferromagnetic exchange 
Ji > and subject to a time-dependent external mag- 
netic field H{t) 
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FIG. 1: T-M /H dc versus T for an Ising chain with Jijk B = 80 
K, in a static magnetic field _Hd c = 1 kOe, calculated for se- 
lected values of the concentration c of non-magnetic impu- 
rities (different symbols refer to c=0.05, 0.02, 0.01 on going 
from bottom to top), and for the pure system (c = 0, full 
line). 



N 



« = -£ 



Jiaj(Tj + i + gn B H(t)aj 



(5) 



III. STATIC PROPERTIES 

For pure CoPhOMe, static magnetization mea- 
surements were originally performed in single- crystal 
samples 1 ^ with a static magnetic field i/dc=l kOe ap- 
plied parallel to the chain direction. A strong increase 
in the quantity T ■ Mj H^c is found on decreasing T be- 
low 100 K, with a maximum reached around 25 K and 
a subsequent decrease. Such a behavior, typical of fer- 
romagnetic and ferrimagnetic systems, was also observed 
in the case of poly-crystalline powder samples of pure 32 
and Zn(II)-doped CoPhOMe^ at different static mag- 
netic fields and different concentrations of Zn(II) non- 
magnetic impurities. It is also a common feature of all 
SCMs identified so faxM 

For an open chain of N spins in an applied field, de- 
scribed by Eq. ([5]), the static properties can be calcu- 
lated analytically in terms of the two eigenvalues, Ao and 
Ai, and of the corresponding eigenvectors of the transfer 
matrix4i2i Denoting the concentration of quenched non- 
magnetic impurities by c, the magnetization per spin of 
the doped chain is 



M, 



doped 



M 



N 



(6) 



and the static susceptibility per spin is 

oo 

Xdopod = E ^ ~ C ) W XN > 



(7) 



JV=1 



where the explicit expression for Mn and \n are reported 
in Appendix A. 

For the infinite chain (i.e. c = and N — > oo), 
the static properties can be calculated analyticall 
terms only of the larger eigenvalue Aq 



3.34 



A = (cosh ft, + \l' 



sinh hn 



-4K 



(8) 



where K = and h = 9 ^^r" ■ The magnetization 
per spin is 



M = k R T 



1 3A 



sinh ho 



-4K 



A clffdc V sinh 2 ho + e~ 
and the susceptibility per spin is 

1 / <9A \ 2 1 a 2 An 



' ( (/ 



X 



k R T 



X 2 \dH d 



(9) 



(10) 



N=l 



In Fig.[TJthe calculated quantity T-M/H^c, see Eq. ©, 
is reported versus T for different values of c (including 
c = 0) and for a fixed value of the dc field (i/dc = 1 
kOe). In Fig. [2 the same quantity is reported for dif- 
ferent values of the static applied magnetic field and for 
a fixed, rather small, value of the non-magnetic impuri- 
ties concentration (c = 0.01). It is worth observing that, 
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FIG. 2: (Color online) T ■ M/H dc versus T for a doped Ising 
chain with Jijks = 80 K, calculated for different values of 
the static applied field, at a fixed value (c = 0.01) of the non- 
magnetic impurities concentration. Inset: experimental data 
for a nominally pure (c = 0) powder sample. 



displayed by the experimental x'i^) f° r T < 10 K, can 
be well reproduced (blue circles) by taking into account 
dynamic effects, at it will be shown in Scction lTVl In sum- 
mary, the new data and calculations reported in Fig. [3] 
not only are in good agreement with the reported pres- 
ence of a distribution of impurities inside the crystal, 26 
but they are also revealing of a non-homogeneous nature 
of the samples. 

A similar two-peaked feature was observed, later, in 
the real part of the ac susceptibility of a different cobalt- 
organic single-chain magnet, for iJdc = 0.5 kOe^i As 
in the case of CoPhOMe, the higher-temperature peak, 
due to the infinite chain, did not present any dynamic 
effect, while the lower temperature shoulder, related to 
finite-size effects, was found to be frequency-dependent. 
In the following we will show that the latter feature can 
be well accounted for by a calculation of the dynamic 
susceptibility in the framework of the simplified model 



at low fields and low temperatures, the curves calculated 
for c = 0.01 resemble the experimental ones, obtained by 
Lascialfari et aL— for a powder sample of nominally pure 
CoPhOMe (see inset). In contrast, according to Eq. ©, 
the theoretical curves for the pure system (not shown 
here) are found to undergo a much stronger increase at 
such low fields and temperatures (as it can also be in- 
ferred from the concentration dependence of T ■ M/H^c 
shown in Fig. Q}. Previous suggestions 2 ^— about the 
presence of lattice imperfections, or impurities that limit 
the chain size even in pure CoPhOMe, appear thus to be 
confirmed. 

In Fig. |3] we show new experimental data for the tem- 
perature dependence of x?(to), the real part of the ac sus- 
ceptibility of a nominally pure CoPhOMe sample, mea- 
sured at frequency v = 1 kHz for a fixed value of the 
static applied field, 7?dc = 1 kOe. For the same field, the 
calculated static susceptibility, x> 01 the infinite chain 
(c = 0, full line) is exponentially vanishing at low tem- 
peratures, see Eq. (|10p. in strong contrast with the di- 
verging behavior of x m zer0 field. For the doped chain 
(c = 0.047, dashed line), the static susceptibility was 
calculated using Eq. (J7J. It turns out that Xdoped, cal- 
culated within the simplified model (JSJ), is in very good 
agreement with previous results from an exact transfer 
matrix calculation^ performed for the more complete 
model Qj as well as with ac susceptibility data in zinc- 
doped CoPhOMe^ 6 - In particular the feature of a peak 
gradually shifting to lower temperatures with increas- 
ing c, a clear signature of finite-size effects, is recovered 
(not shown here) by our simplified model. From Fig. [3] 
one sees that the two-peaked feature, experimentally ob- 
served at nonzero field in x'( w ) versus T, can be well 
reproduced by taking a weighted average (light-grey tri- 
angles) between the static susceptibility of the pure chain 
(black full line) and that of the doped one (black dashed 
line). Finally we mention that also the steep decrease, 



IV. DYNAMIC PROPERTIES 

In ac magnetic measurements, a small ac drive mag- 
netic field is superimposed on the dc field, causing a time- 
dependent moment in the sample. Therefore we are faced 
with the theoretical problem of determining how the ki- 
netic equations of motion for the time-dependent spin av- 
erages of a finite, open chain in zero magnetic fiel d 21 ' 22 
are modified by the presence of magnetic field of the gen- 
eral form 

H(t)=H dc + H 1 e^ t (11) 

i.e. the sum of a static dc field of any intensity, H dc , 
and of an ac field, oscillating at the angular frequency u. 
In the following we will make use of the reduced fields 

ho = 9 ^"t c and h i = 3 k%T k - As the experimental 
oscillating fields are usually much smaller than the static 
field, we will consider the case Hi <C -Hdc, which allows 
us to use the expansion (h(t) = gM ^^ t - > ) 

tanh/i(t) pa tanh/i + hi e~ lut (l - tanh 2 ho). (12) 

Generally speaking, the susceptibility induced by a field 
as in Eq. fjl If) will have a real and an imaginary part 

xM = x» + ix» (13) 

from which the dynamic behavior of molecular materials 
is usually extracted. In the following we will thus focus 
on calculating these experimentally relevant quantities, 
which are compared to data in the next section. 

Typical values of the frequencies, v = used in 
ac susceptibility measurement s] 10 ' 26 on pure and doped 
CoPhOMe, range between 0.1 and 10000 Hz,. In con- 
trast, one has v sa MHz for proton nuclear mag- 
netic resonance (NMR) and muon spin relaxation (/iSR) 
experiments! 35 ' 36 In the following, before passing to the 
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FIG. 3: (Color online) (a) Real part of the ac susceptibility, x'fa): (square symbols) measured versus temperature in a 
nominally pure CoPhOMe sample for frequency v = 1 KHz and fixed static applied field Hdc = 1 kOe. (b) Static susceptibility, 
X, calculated with Ji/ks — 80 K and Hdc = 1 kOe, for a doped Ising chain (c = 0.047, dashed line) and for an infinite chain 
(c = 0, full line). The two-peaked structure (light-grey triangles) is obtained by taking a weighted average of the two quantities 
(80 % versus 20 %, respectively). The real part of the dynamic susceptibility, X ( u )j calculated according to the theory in 
Section HVl is denoted by blue circles. 



comparison with experimental data, we first examine the 
properties of an infinite chain (subsection A), then in- 
clude finite-size effects (subsection B), and eventually ex- 
plore the single-frequency approximation (subsection C). 



A. Infinite chain 

Before considering finite-size effects, it is instructive 
to calculate first the ac susceptibility of an infinite Ising 
chain in presence of a static magnetic field. The kinetic 
equation for the site-independent average m(t) — (o~(t)} 
of a spin in the infinite chain is 



ro^p = -(l-7)m(i) + 



l-iTi(t) tanh/j(t) (14) 



where 7 = t&nh(2K) and Ti(t) is the nearest neighbor 
spin-pair time-dependent correlation function. The ki- 
netic equation for Ti(t), on its turn, involves higher order 
time-dependent correlation functions, so that eventually 
an infinite sequence of equations is obtained as a conse- 
quence of h(t) ^ 0. In order to truncate this hierarchy, 
we adopt the local-equilibrium approximation ! 29 ! 30 i.e., 
the relation holding at equilibrium^ between the near- 
est neighbor correlation function and the magnetization 
(where m eq is defined in Eq. |9|), 



"-eq 

(1 - Kq) 



2 cosh hg — \J sinh 2 h + e 4K 
cosh ho + v s i n h + e~ 4K 



(15) 



is assumed to hold locally also at any time t ^ 0. The 
main advantages of this approximation, first proposed by 
Huang^S for an infinite Ising chain model with Glauber 
dynamics in H^ c =/= 0, are: (i) it provides the exact 
steady-state solution^ - in contrast with the mean field 



approximation, which assumes simply T% , 



e 9 ; (ii) 



it is valid for any value of the applied field, in contrast 
with the perturbation method, which assumes ho <g; l. 2 
In this way, a nonlinear equation for m(t) is obtained, 
where the approximation of a linear response of the 
chain applies only to the ac field, H\. We thus assume 
6m(t) = m(t) — m eq , i.e. small departures of the mag- 
netization from its equilibrium value m eq . Likewise, we 



expand T±(t) » Ti\eq 
account that ^$\ eq 



^$j\eq Sm(t) and, taking into 
2 tanh ho at equilibrium^ we fi- 



nally obtain a linear non-homogeneous differential equa- 
tion for Sm(t) 



d Sm(t) , , 9 , , _ , „ 

T Q = -(1 - 7 + 27tanh 2 ho) Sm(t) 



dt 
+hi e 



1 - tanh 2 ft Q )(l-7 I\ 



,eq 



(16) 



In absence of the ac field (hi = 0), one finds an expo- 
nential time decrease for the magnetization fluctuation 



8m(t) = 5m{to)e~ x ^ [t - to)lT °. 



(17) 



Thus, for the infinite Ising chain, there is a single relax- 
ation time, Too, related to the a-dimensional parameter 
Aoo by 

Aoo = — = 1 -7 + 2 7 tanh 2 /i . (18) 

Too 

Notice that, for ho — > 0, Glauber's result 2 - of an expo- 
nentially diverging relaxation time at low temperatures 
( T oo ~ \ T o^ AK for k B T -c Jj) is correctly recovered, 
while for H^c 7^ the relaxation time of the infinite chain 
does not diverge any more. 

In presence of the ac field (hi ^ 0), the general solution 
of Eq. (fTU) is 



5m(t) = 6m(i )e- A ~ (t ~ to)/ro + hi (1 - tanh 2 ho) 



6 



(l-7ri,eg) / e^'V^'"*'^ dt'. (19) 



in 



Taking into account that Aoo 7^ 0, one can safely let to — > 
— 00 in order to find a solution that does not depend on 
the initial conditions 2 * 3 ^ 



Sm(t) = (1 - tanh 2 h )(l - 7 T 1}eq ) '' ' ' 



(20) 



The fluctuation of the total magnetization of the infinite 
fsing chain is obtained summing over all the N spins and 
letting N — » 00 

5{M(t)) = Ng» B 5{a(t)) = X ^)H ie - iut , (21) 

so that the dynamic susceptibility, xi^i IS 

a 2 u 2 N 

X (cj) = (1 - tanh 2 h )(l - 7 T heq ) 

K B 1 

1 Xoo 



Aoo — ilotq 1 — iuj-p- 



(22) 



where Xoo denotes the static susceptibility of the infinite 
Ising chain for H^c 7^ 

= 9 2 ^ B N (1-tanh 2 fe )(l - 7 I\ eg ) 
X °° k B T 1-7 + 2 7 tanh 2 h ' 

For /i -> 0, Glauber's result 

g 2 »%N I + 77 1 



(24) 



is correctly recovered, by taking into account that A~ 

2rj_ 



1 — 7, ri, e9 — > r\ = tanh if, and 7 = jr 2- 



B. Finite chain 

In the case of an open Ising chain with a finite num- 
ber N of spins, the lack of translational invariance leads 
to N kinetic equations for the N site-dependent spin av- 
erages 5{a p (t)} (p — 1, • • • , N)£^£2£2- As in the case of 
the infinite chain, we can introduce the local-equilibrium 
approximatio n 29 ' 30 in order to truncate the infinite se- 
quence of equations for the higher-order time-dependent 
spin correlation function. Next we perform the lineariza- 
tion of the kinetic equations, in the hypothesis of a linear 
response to the ac magnetic field. We then obtain a set 
of N linear differential equations, which can be written 
in matrix form 

t — = -Y ■ £ + h ie - iujt (l - tanh 2 ho)®. (25) 
dt 

where £ and '3/ are N x 1 vectors containing the spin fluc- 
tuations and the non-homogeneous terms, respectively 
(see Appendix B for details). Y is a real, symmetric, 
tridiagonal, N x N matrix, with nonzero a-dimensional 



eigenvalues Xj (j — 1,...,N), while <I>( a j) are the corre- 
sponding N X 1 eigenvectors. In the limiting case ho — > 0, 
the numerical solutions for Xj coincide with the ones 
obtaine d 21 ' 22 in the framework of a finite-size scaling cal- 
culation of the Glauber dynamics in a zero static field. In 
particular, the low-temperature expansion (k B T <C Ji) 
for the eigenvalue of a finite open chain with N spins 

is Xi(Hdc = 0) « -^2_e k B T j2i to be compared with 
-Hi 

Aoo(iidc = 0) « 2e k s T for the eigenvalue of the infinite 
chain <2 

We are interested in the long time behavior of the 
system, characterized for being independent on the ini- 
tial condition. Thus, solving Eq. (I2"5j) by the method of 
eigenfunctions 3 - 9 . and letting to — )• —00^^ we obtain for 
the fluctuation of a single-spin average (p = 1, • • • , N) 

6(a p (t)) = hie- iut (l-ta,uh 2 ho) 



N 



IUJTq 



N 

^(E^*™) ( 26 ) 



The fluctuation of the magnetization of the finite open 
Ising chain is obtained summing over the N spins 

JV 

S(M N (t)) = 9 fi B *<<>■*(*)> = XN(u>)H ie - iut (27) 
P =i 

The dynamic susceptibility xn(u) takes the form 

xn{uj) = n 2 + J Mx ^ HAc) (28) 

where the angular frequencies are 

A; 



T 



(29) 



and the corresponding frequency weights are 

2 AT 
1 — tanh h 



Aj(Xj,T,Hdc) 



A, 



P =i 



A' 



E$( A i)lIl 
m r 



m— 1 



(30) 



The angular frequencies (I2U1) are expressed in terms of 
the a-dimensional eigenvalues, Aj, of the matrix Y and of 
the characteristic time, To, for the spin flip of an isolated 
spin. The latter is a free parameter which is expected 
to depend, in general, on the intensity of the applied 
magnetic fields The general expression for the dynamic 
susceptibility per spin of a doped chain is thus 



N=l 



(31) 
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H =2 kOe N=11 

do 
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FIG. 4: (Color online) Temperature dependence of the first three eigenvalues Aj = ttjTo (j = 1, red full circles; j = 2, 
green full squares; j = 3, blue full triangles) of an open Ising chain, calculated for a finite number of spins (iV=ll and 
N = 101) in presence of a nonzero static field (Hdc = 2 kOe and Hdc = 10 kOe). For comparison, the eigenvalue Aoo 
(black open circles) of an infinite chain in nonzero field is also reported. Dashed lines denote the zero-field, low-temperature 

expansions for the smallest eigenvalue of a finite chain, Ai (Hdc = 0) « N 2 _ 1 e k B T ,— and for the eigenvalue of an infinite chain, 

Aoo (Hdc = 0) ~ 2e fc s T .— The horizontal lines, from bottom to top, denote the quantity 2-kvtq, calculated for three different 
frequencies: v = 1 kHz, ^(MHz)= 4.26 ■ i?dc(kOe), and v(MHz)= 13.55 ■ //dc(kOe), typically used in ac susceptibility, proton 
NMR, and /xSR measurements, respectively. Insets: calculated temperature dependence of the frequency weights, Aj(Xj,T, Hd c ) 
with j = 1 and j = 3 (red full circles and blue full triangles, respectively); for odd N, the weights of the even modes are zero^S. 
The curves denoted by black open circles and black dashed lines represent the frequency weights of an infinite chain in fide 7^ 
and in Hdc = 0, respectively. The parameters used for the calculations were Ji/Ub = 80 K, g = 2, and to = 4 ■ 10 -13 s. 



C. Single-frequency approximation 



frequency (uj = 0) in Eq. 



N , rp 

^A ] (\ 1 ,T,H dc )= 11 



3=1 



XN 



(32) 



Let us start deriving a quite general sum rule for the where \n, the static susceptibility of a finite open Ising 
frequency weights, which is readily obtained letting zero chain with N spins, subject to the dc field H& c , is given 
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FIG. 5: (Color online) Temperature dependence of the calculated dynamic susceptibility for an infinite (c = 0) ferromagnetic 
Ising chain in zero static field Hdc = (figures a and b), and of a doped (c = 0.047) Ising chain in nonzero static field Hdc = 2 
kOe (figures c and d). The different curves refer to different frequencies (u — 0.1, 1, and 10 kHz on going from left to right) of 
the ac magnetic field. The calculated static susceptibility of the infinite chain in Hdc = 0, is denoted by the full black line: it 
diverges exponentially as T is decreased to 0, see (a), while in Hdc 7^ it goes through a peak, and then vanishes exponentially, 
see (c). The calculated static susceptibility of the doped (c = 0.047) chain in nonzero (iidc = 2 kOe) field is denoted by the 
black dashed line, see (c). 



by Eq. (|A5|) . Next we observe that, in the approximation 
of a single characteristic frequency, fl c , dominating the 
relaxation of the magnetization fluctuations, the equa- 
tion (|28p for the dynamic susceptibility of a finite, open, 
Ising chain with ./V spins assumes the simple form 



Xjv(w) ~ XN 



fl 2 + ioj£l c 

fl 2 +UJ 2 



(33) 



In principle, the characteristic frequency is not necessar- 
ily fl c = Xi/tq, i.e. related to the smallest eigenvalue 
of Y . Rather, in order to determine the dominating 
frequency, the temperature dependence of the frequency 
weights must be taken into account. 

In the single-frequency approximation, the dynamic 
susceptibility per spin of a doped chain also simplifies 
greatly 



XM ~ Xdopcd 



Q 2 + iujflc 

n 2 + oj 2 



(34) 



where Xdopod is given by Eq. 0. 

Finally we observe that, using the fluctuation- 
dissipation theorem, the linear response Sn(uj) of a fi- 
nite open Ising chain with N spins can be expressed as a 



weighted sum of N Lorentzian functions centered at zero 
frequency with widths ttj 



S n (uj) 



2k B T „ 

N 

2 sVIE cfi±7a M^T,H Ac ). (35) 
j=i 3 



When the approximation of a single dominating fre- 
quency (|33|) holds, the linear response is simply the prod- 
uct of T\n and of a single Lorentzian function, centered 
at zero frequency, with width equal to the characteristic 
frequency VL C 



S N (u>) = — xM » 2k B T XN 



a 



<> 2 



(36) 



40 



It is worth noticing that expressions quite similar to 
Eqs. (|35|) and (|36|) were obtained for the spectrum of flue 
tuations of a cluster magnetization by Santini et 
in the framework of an exact calculation of the energy 
levels of three important classes of magnetic molecules 
in contact with a phonon heat bath: namely, antifer- 
romagnetic rings, grids, and nanomagnets. Moreover, 
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FIG. 6: (Color online) Temperature dependence of the measured ac susceptibility of CoPhOMe. The data refer to the same 
sample as in Fig. [3ja). In (a) and (b), the frequency was fixed (y = 1 kHz) and the dc field was varied (Hdc = 0.333, 0.666, 
and 1 kOe on going from top to bottom). In (c) and (d), the dc field was fixed (-Hdc = 1 kOe) and the frequency was varied 
(y = 10, 100, and 1000 Hz on going from left to right). 



Bianchi et al^ recently showed that, while for antifer- 
romagnetic homometallic rings the approximation of a 
single dominating frequency is valid, it does not hold for 
heterometallic rings^ due to the presence of inequiva- 
lent ions which prevent mapping local-spin correlations 
with the corresponding total-spin ones. 



spin flip of an isolated spin, r , was left as a free pa- 
rameter (see later on). Moreover the field dependence of 
ro, though expected in principle^ was neglected for the 
sake of simplicity. 



V. RESULTS 

In this section we compare the theoretical results de- 
rived in [IV] with experimental ac susceptibility data, ob- 
tained in nominally pure and zinc-doped CoPhOMe for 
different values of Hdc and of the frequency of the oscil- 
lating field. Also, new experimental data for the relax- 
ation time of the magnetization, measured as a function 
of the static magnetic field at fixed temperature, will be 
discussed and compared with theoretical calculations. 

Let us first provide evidence for the correctness of the 
single frequency approximation, Eq. (|34j) . by showing the 
temperature dependence of the eigenvalues and weights 
of a finite, open, Ising chain in an applied dc field. The 
value of the exchange constant we assumed for the cal- 
culations, J 1 /kg — 80 K, is knownSSi to provide the cor- 
rect temperature dependence for the static and dynamic 
properties of CoPhOMe. The characteristic time for the 



In Fig.@]some of the N calculated a-dimensional eigen- 
values, Xj, of the real tridiagonal matrix Y, defined in 
Eq. (|25|) and Eq. (|B10[) . are reported as a function of 
inverse temperature for different values of the number 
of spins, N, and of the applied static field, H^. For 
the sake of comparison, also the temperature dependence 
of the eigenvalue Aoo, see Eq. (fl~8)l . of an infinite chain 
in the same field is shown. Except for the case of long 
chains and strong fields, at sufficiently low temperatures 
a single mode dominates the low-frequency dynamics of 
a finite, open, Ising chain with N spins: namely, the 
mode with characteristic frequency f2 c = Xi/tq, where 
Ai is the smallest eigenvalue of Y. The temperature de- 
pendence of the frequency weights corresponding to the 
various modes is displayed in the insets of Fig. [4] for not 
too long chains and not too strong fields, one can see that 
the frequency related to the smallest eigenvalue Ai has 
the strongest weight. 
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FIG. 7: (Color online) Temperature dependence of the dynamic susceptibility of CoPhOMe, calculated for a doped chain 
(c = 0.047) in the single-frequency approximation, see Eq. (|34p and Eq. lff|l. for the same values of frequency and dc field as in 
Fig. [6] Thin dashed lines denote the static susceptibility of the doped system (c = 0.047) versus T. Thin solid lines denote the 
dynamic susceptibility of the infinite chain (c = 0): notice the absence of any frequency dependence for fixed H,ic = 1 kOe, see 
(c), at the chosen values of v. 



A. Ac susceptibility 

Some years ago, the ac susceptibility of CoPhOMe was 
measured versus T in single crystals, both for H^c = Oe 
in a nominally pure sample 1 - and for H^c = 2 kOe in a 
doped one. 26 In this work, we present some new data for 
x(u;) versus T. Data were obtained, using an homemade 
ac probe and a Cryogenics magnetometer, on a collec- 
tion of nominally pure single crystals of CoPhOMe. The 
crystals were all aligned with the chain axis along the 
magnetic field direction. The frequencies in an ac sus- 
ceptibility experiment typically range between 0.1 and 
10 kHz, and we can safely adopt the single-frequency 
approximation in order to account for the temperature 
dependence of m moderate fields (-f/dc < 2 kOe). 
Additionally the oscillating fields used were always be- 
low 8 Oe, and the condition H dc 3> H\ is also fulfilled. 

In Fig. [5fa)(b) we show the temperature dependence 
of the real and imaginary part of the dynamic suscepti- 
bility calculated for an infinite chain with H& c = 0, while 
analogous quantities calculated for a zinc-doped chain 
(c = 0.047) in H dc = 2 kOe are reported in Fig. EJc) (d) . 

In the pure system (c = 0) with H^c — 0, see Eq. (|2~4"1) . 
a single, resonating peak is found both for x'(w) and 
x"(w) versus T, whose position gradually shifts to higher 



temperature with increasing the frequency of the tiny os- 
cillating field, see Fig. [5ja) (b) . The phenomenon can be 
interprete d 3 ^ 42 ' 43 as a manifestation of stochastic reso- 
nance in a set of coupled bistable systems: i.e. there is 
an optimal value of noise, for which the response of the 
dynamic system to the driving field is maximum. In a fer- 
romagnetic or ferrimagnetic chain, the role of stochastic 
noise is played by thermal fluctuations, and a resonance 
peak occurs when the deterministic time scale of the os- 
cillating magnetic field matches the statistical time scale 
associated with the spontaneous decay of the net magne- 
tization. 

For the doped system in H dc = 2 kOe, we find that 
a frequency-dependent peak in the calculated x'( w ) an d 
x"(w) versus T (the colored curves in Fig.[5jc)(d)) devel- 
ops at substantially lower temperatures with respect to 
the peak in the static susceptibility of an infinite chain 
(the full black curve in Fig.JSfc), see Eq. E5|) . This can be 
easily understood by looking at Fig. [4] for a finite, open 
Ising chain in a moderate field, the fulfillment of the res- 
onance condition (w = J7 C in Eq. I34p . occurs at low tem- 
peratures, as signaled by the crossing between the full 
horizontal line (which represents a typical value, 1 kHz, 
of the frequency v in an ac susceptibility measurement) 
and the curve (red full circles) representing the T depen- 
dence of the smallest frequency mode, X\. In contrast, 
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FIG. 8: (Color online) (a) Frequency dependence of the imaginary part of the ac susceptibility, \" , measured in CoPhOMe at 
fixed temperature T — 9 K for selected values of the static magnetic field, (b) Field dependence of the relaxation time r of the 
magnetization of CoPhOMe as deduced from the peak position of x" versus frequency at T = 9 K. Squares denote experimental 
results; the lines connect theoretical values calculated for a doped chain in the approximation of a single dominating frequency 
(see Eq.|M]and Eq.0, using Ji/k B = 80 K, g = 2, and t = 4 • 10 -13 s. 
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FIG. 9: (Color online) (a) Time dependence of the magnetization m(t), measured in CoPhOMe at fixed temperature T = 5.1 K 
for selected values of the static magnetic field, (b) Field dependence of the relaxation time r of the magnetization of CoPhOMe. 
Squares denote experimental results obtained by an exponential fit of m(t); the line connects theoretical values determined as 
in Fig. 7b, but for T = 5.1 K. 



the crossing does not occur in the case of an infinite chain 
(black open circles) subject to the same dc field: i.e. the 
relaxation rate of the infinite chain does not fulfill the 
resonance condition (ujtq — Aoo in Eq. I22|) at the low 
frequencies involved in the ac susceptibility experiment: 
thus, the infinite chain does not present any dynamic re- 
sponse. As already observed in Sect. Ill, nominally pure 
samples are non-homogeneous and consequently, at low 
temperatures, only the regions with dilute chains con- 
tribute to the dynamic properties in a significant way. 
This explains the experimental result o 10 ' 11 ! 23 " — in which 
the measured relaxation rate, for H=0, was always found 
to follow the modified Arrhenius law with a halved en- 



ergy barrier, r(T) ~ tq(L) e fc s T , where L is the mean 
length of a chain. 

In Fig. [B] we present new experimental data for 
the temperature dependence of the ac susceptibility of 
CoPhOMe, measured for a fixed value of the frequency 
(v = 1 kHz) at different values of the static magnetic 
field (-ffdc < 1 kOe), see (a)(b), and for a fixed value 
of the dc field (Hdc = 1 kOe) at different values of the 
frequency (v < 1 kHz). In Fig. [7] we plot the real and 
imaginary part of the dynamic susceptibility, calculated 
using Eqs. ((Ml) and © for a doped chain (c = 0.047) at 
just the same values of frequency and dc field as in Fig. 
For the sake of comparison, also the dynamic susceptibil- 



12 



ity of the infinite chain (c = 0) was calculated: for fixed 
-ffdc = 1 kOe, at the chosen values of v, this quantity 
was not found to display any frequency dependence (see 
the black line in Fig.[7][c)), in fine agreement with exper- 
imental observations (see the higher temperature peak in 
Fig.EJc)). Thus we can conclude that, for the moderate 
values of v and H& c exploited in the present ac suscepti- 
bility experiments: (i) the frequency dependence of x( w ) 
can be attributed solely to finite-size effects; (ii) the ap- 
proximation of a single dominating frequency is a good 
one for the calculation of the dynamic susceptibility. 



B. Relaxation time 

In Fig. Elja) we show the frequency dependence of the 
imaginary part of the ac susceptibility, measured for a 
sample of five crystals of nominally pure CoPhOMe, at 
a fixed temperature of 9 K and for iJdc ranging from 0.1 
to 2 kOe. Fitting of the curves was performed using an 
extended Debye model to extract the peak frequency^ 



xM = Xs + 



XT - XS 
1 + (iujt) 1 - 



(37) 



where xt is the isothermal susceptibility, \s is the adia- 
batic susceptibility, and a is a parameter that measures 
the distribution of relaxation times in the sample. In the 
present measurements it was found that a ss 0, meaning 
that the whole system relaxes with a single characteristic 
time. The so-obtained field dependence of the relaxation 
time, t, is reported (red squares) in Fig.[SJb). 

The field dependence of the relaxation time at a lower 
temperature, T = 5.1 K, was instead obtained perform- 
ing dc measurements. The magnetization was first sat- 
urated in = 25 kOe, the field was swept to a new 
Hdc value with a field sweep rate of about 250 Oe/s, and 
the time dependence of the magnetization of the system 
was recorded in presence of an applied static magnetic 
field. The relaxation of the magnetization as a func- 
tion of time in nonzero field, measured using a Cryo- 
genic S600 superconducting quantum interference device 
(SQUID) magnetometer, is reported in Fig. Ufa). The 
magnetization was fitted using a mono-exponential law, 
M(t) = M(t )+at e ~(*- to )/ T , and we obtained the relax- 
ation time t versus H reported (red squares) in Figure 

For the interpretation of the r versus H data, con- 
centrations of non-magnetic impurities in the range c ss 
0.02 — 0.05 were assumed, i.e. comparable with those 
evidenced by ac and dc susceptibility measurements in 
the same compound. The theoretical curves for r versus 
H, reported as full (c = 0.02) and dashed (c = 0.047) 
lines in Figure EJb) and EJb), were determined as fol- 
lows. The dynamic susceptibility x( w ) was calculated 
in the approximation of a single dominating frequency, 
see Eq. (|33[) . and the frequency maximum of x'{ u ) was 
taken to determine r = 1/ (w max ) . In order to reproduce 
the correct order of magnitude for the relaxation time r 



of the magnetization of CoPhOMe (both at T = 5.1 K 
and 9 K), one has to assume a value r « 4 • 10~ 13 s 
for the characteristic time of spin flip of an isolated spin 
(the only free parameter in Glauber's theory^). For the 
sake of simplicity, the same field-independent value was 
assumed in all calculations throughout the paper. 

As regards the dc field dependence of the relaxation 
time of the magnetization, t, we can conclude that it is 
well reproduced by theory, provided that finite-size ef- 
fects are included. In fact, in a measurement performed 
for an inhomogeneous CoPhOMe sample, no apprecia- 
ble contribution to the relaxation time, r, is expected 
from the infinite chain at low temperatures, owing to the 
high value of the exchange coupling (Ji/ks = 80 K). 
For example from Eq. (Tl5|) . putting ro = 4 ■ 10 -13 s, at 
T = 5 K one has ToD = 1.25 • 10 15 s for H dc = 0, while 
in a static field H dc — 10 Oe the relaxation time reduces 
to Too = 2.77 • 10~ 6 s. In contrast, an experimentally 
appreciable relaxation time is associated with finite-size 
chains. 



VI. CONCLUSIONS 

In this paper, the spin dynamics of the archetypal 
molecular magnetic chain, CoPhOMe, in presence of an 
external magnetic field of any intensity, was investigated 
using a simplified model, consisting of a one-dimensional 
Ising ferromagnet with a stochastic dynamics caused by 
the interaction of the spins with a heat reservoir^ 

In the framework of a local-equilibrium approximation, 
devised to truncate the infinite sequence of kinetic equa- 
tions originated by the presence of a nonzero static mag- 
netic field, and of a linear response of the system to a 
small oscillating field, we first calculated the dynamic 
susceptibility of an infinite chain. Next, the theory was 
generalized to a finite, open Ising chain. We showed that 
the dynamic susceptibility of an open chain with a fi- 
nite number TV of spins can be expressed as a weighted 
sum of N frequency contributions, related to the N re- 
laxation rates of the magnetization fluctuations. From 
the comparison with the ac susceptibility data obtained 
for nominally pure samples we can draw two conclusions: 
i) the pure samples are really non- homogeneous, because 
regions with very low density of defects coexist with re- 
gions with relevant density of defects; ii) only the latter 
regions show a contribution to the dynamic relaxation, 
in the frequency range conventionally investigated by ac 
susceptibility measurements. 

For doped CoPhOMe chains, the approximation of a 
single dominating frequency was found to be quite sat- 
isfactory in order to account for ac susceptibility data— 
in a moderate static magnetic field (-ffdc = 2 kOe), both 
as a function of temperature and of the frequency of the 
small oscillating field. 

Finally, it is worth observing that, on the basis of our 
calculation of the T- and ^-dependence of the relaxation 
rates (and corresponding frequency weights), we do not 
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expect that the approximation of a single dominating fre- 
quency will be able to account for X H nuclear magnetic 
resonance (NMR) and muon spin rotation (/iSR) experi- 
mental data22i22i2§j2£ in pure and zinc-doped CoPhOMe, 
because these techniques probe the local spin dynamics 
at frequencies (v ~ MHz) substantially higher than the 
typical frequencies (y ~ kHz) used in an ac susceptibil- 
ity experiment. Anyway, the present results constitute a 
fundamental background for the future theoretical inves- 
tigation of such high-frequency regimes. 
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Appendix A: Analytical expressions of the magnetization AIn and the static susceptibility \n of an open 

Ising chain of TV spins in a dc magnetic field Hdc 

Following Wortis,— the free energy Fn of a finite chain of N spins, coupled by the Hamiltonian ([5]) with H = H^c 
and subject to open boundary conditions, can be expressed as the sum of a bulk, a surface, and a finite-size contribution 



Fj 



N 



k B T{ AlnA + ln 
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In 
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L (-) 
o VAoJ 



(Al) 
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Ao, 



3 K ( cosh h ± \/sinh 2 h + e~ iK 



(A2) 



Aq i = cosh h ± 



sinh 2 h + e~ 2K 
Vsinh 2 h + e~ iK 



(A3) 



with K = jf^f, ho = g ^y dc ■ For the reader's convenience, explicit expressions for the static magnetization, Mjv, 
and static susceptibility, xn, of a finite open Ising chain with N spins in presence of a magnetic field are reported 
hereafter. The magnetization is 
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(A4) 



and the static susceptibility is \n = — 
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(A5) 



The derivatives with respect to the magnetic field in Eq. (IA4[) and Eq. (| A5|) are expressed as 
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Appendix B: Relaxation times of an open chain of N Ising spins in an applied magnetic field 



In order to calculate the relaxation times of an open chain of N spins coupled by the Ising exchange Hamiltonian 
and subject to a static magnetic field H, Eq. ([5]), following Glauber— one starts by writing the kinetic equation of 
motion for the time-dependent average of an interior spin ((<7 p ) with 1 < p < N) 



d(<Jp) 
dt 



+ o ( (°p-i) + ( ct p+i) ) + tanh/i(i) 1 - ~ ( {<j p -\cr p ) + (a p <r p+1 ) 



(Bl) 



where To is the characteristic time for the spin flip of an isolated spin, 7 = tanh(2i^T) and hit) = ■ For the two 

spins (p = 1 and p = N) located at the ends of the open chain, following Coulon et aZ.— one has instead 



to- 



d(ai) 
dt 



= — (ci) + 77(02) + tanh/i(t) 1 — 77(01 02) 



(B2) 



to dt = -(cat) + Vi^N- 



}+tanh/i(t) 1 — rj{<TN-i<JN) 



(B3) 



where r\ = tanh K£ In our case, see Eq. (jTTJ) , the time-dependent magnetic field is assumed to be H(t) = H^ c +Hi e lut 
i.e. the sum of a static dc field of any intensity, i/do an d of a much lower AC field, Hi <C -ffdc, oscillating at the 
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angular frequency ui. Thus we can expand tanh h(t) ss tanh/io + iut (l — tanh 2 /io), where ho = 9M fc B y dc and 
111 ~ fe B T • 

Following Coulon et al.^ a linearization of the equations (|B1[) . (|B2|) . and (|B3I) around the spin equilibrium values 
is next performed in the framework of linear response theory (i.e. small departures from thermal equilibrium are 
assumed) 



(a p ) « {<Tp)N,eq + 5(a p ) 

(a p a p+ l) S» {(Tp(Jp+l) N ,eq + S((7 p a p+ i) 



(B4) 



In this way, we obtain a system of N differential equations in the spin fluctuations (S(a p ) with p = 1, • • • . AT) where, on 
the right hand sides, the presence of the field involves the presence of variations of time-dependent nearest neighbors 
spin-spin correlation functions {5{a p a p +i)) 



To- 



TO 



dd(ai) 

dt 
dS(a p ) 



df 

d 6(cr N ) 
dt 



= — 5((Ti} + rj 8{u2) — tanh/io S(a 1 a 2 ) + h x e %U3t (\ - tanh 2 h ) l~rj(a 1 a 2 )N. 
= -5{<Tp) + ^ (6(ap-i) + 5{a p+ i)^j - t&nh h {a p-xap) + 5(a p a p+ i)^j 

— ' t 2 r t ( 

+ hie tul (1 - tanh ho) 1 - - 1 (<r P -i<J P )N,eq + (cr p a p+ i) eq 



= — (5(<7Ar)+?7 6(<tn-i) — tanh ho S^n-i^n) + hie tu>t {\ — tanh 2 ho) 1 — r/^N-i^N) N,eq (B5) 



If one writes down the kinetic equations for 8{a p a p+ i) , one finds and infinite sequence of equations involving other, 
higher-order spin correlation functions. Such an infinite hierarchy of equations can be decoupled by resorting to the 
local-equilibrium approximation first proposed by Huang in the case of the kinetic equation of an infinite Ising chain 
in a magnetic field:' 5 ' i.e. the relation existing between the magnetization and the correlation function at thermal 
equilibrium^ is assumed to hold locally although the system is not in equilibrium (i ^ 0). The advantage of the local- 
equilibrium approximation with respect to perturbative methods^ or mean field approximation is that it provides an 
exact steady-state solution.— 

In the case of a finite chain with N spins, the local-equilibrium approximation was implemented, following Coulon 
et al.j22 by expressing the variations of two-spin correlation functions in terms of a linear combination of the variations 
of single-spin averages 



The coefficients are 2 - 

An,p — 



5{a p a p+ i) = A NtP 6{o- p ) + B NtP 5(o- p+ i) 

(°~N-p) N-p,eq ~ V (< J p)p,eq 



q _ (°~p)p,eq ~~ V (o'N-p) N-p,eq 
1 — V (o~p)p,eq{&N-p)N-p,i 



1—1] (<Tp)p je g(o'JV— p) N—p,eq 1 7 "," ;>..</',". \ 

where the average spin values can be calculated, at thermal equilibrium, by using the recursive relation 45 

tanh/io + r\ {a p ^x) p -^ eq 



eq 



l + i] tanhft-o (op-i)p-i, 



(B6) 



(B7) 



(B8) 



with the initial condition {a\)i , e<3 = tanh/i . The set of iV linear differential equations in the N variables 5{a p ) can 
be rewritten in matrix form as 



hie~ iult {l - tanh 2 h )V. 



(B9) 



The x 1 vector £ has elements S p = S(a p ) (p — 1, N). Y is a real, symmetric, tridiagonal, N x N matrix whose 
nonzero elements are (1 < p < N) 



Yi,i = 1 + V An,i tanh/io, Y± i2 = V (^Bn,i tanh/i -l), 



y p-i,p — 2 \ A N,p-i tanh/iQ-lJ, Y P:P = 1 + - [A NiP + Sjv.p-iJ tanh/i , F p 

Yn,n-i = V \ A NiN -i tanh/iQ-lJ, Fjv.jv = 1 
The JVxl vector * in Eq. (|B9|) has elements 



+i = ~ (-B„ iP tanh h - 1 
?7 Sjv,jv-i tanh /i 



(BIO) 



1 -ry (<Ticr 2 )iv,. 
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\Pp-lPp)N,eq + (OpCTp+l) N ,eq 



1 < p < N 



1 - V {^N-l(7 N )N,eq = *1- 



(Bll) 



This equation can be solved, e.g. using the method of cigcn functions, 39 as described in Section IV. B. 
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